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ABSTRACT 


This  thesis  studied  the  impact  damage  tolerance  of 
composite  cylindrical  panels  subjected  to  low  velocity 
impacts.  Graphite/Epoxy  and  Kevlar/Polyester  panels, with 
ply  lay-ups  of  [0/90]o  .  [±45/90/01,  and  [0/90]-  >w ere 
(^investigated. 


STAGSC-1 ,  the  finite  element  program  used  to  perform 
the  analysis,  can  not  analyze  a  laminated  woven  material 
such  as  Kevlar/Polyester,  so'S^simple  method  was  devised  to 
model  the  material  as  an  orthotropic  laminate.  The 
equivalent  material  properties  agreed  to  within  5*  of  the 
experimental  properties  determined  by  standard  material 
property  tensile  tests. 

%  nonlinear  collapse  analysis  was  performed  to  predict 
the  strength  of  the  undamaged  panels.  The  experimental 
collapse  loads  differed  by  as  much  as  62*  (lower)  from  the 
finite  element  predicted  loads.  ""The  sensitivity  of  the 
composite  panel  to  small  imperfections  accounts  for  the 
difference  between  analytical  and  experimental  loads. 
Imposing  a  transverse  sinusoidal  imperfection  with  a  maximum 
amplitude  equal  to  the  thickness  of  one  ply  resulted  in 
analytical  loads  that  agreed  with  experimental  loads  to 
within  16*  at  worst. 


Graphite/Epoxy  and  Kevl ar/Pol yester  curved  panels  were 
impacted  with  energies  of  1-3  ft-lbs  (velocities  of 
3.40-5.36  ft/sec).  No  damage  occurred  for  impact  energies 
below  1-3  ft-lbs,  depending  on  the  ply  lay-up.  The  residual 
strength  of  the  damaged  panel  was  only  2-6*  lower  than  the 
strength  of  undamaged  panels. 

Ultrasonic  C-Scan  and  microphotographs  were  taken  of 
panels  impacted  with  energies  between  8-17  ft-lbs 
(velocities  of  9.62-13.89  ft/sec)  to  investigate  and 
characterize  impact  damage  from  the  material  response  point 
of  view.  Even  at  low  impact  energies,  C-scans  and  stereo 
X-rays  show  that  internal  damage  occurs  even  though  no 
exterior  signs  of  damage  are  visible.  The  damage  can  be 
characterized  by  matrix  cracking,  delamination,  and  fiber 
breakage.  ^  ;>  *  .  •  •  x  •  (  > 
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IMPACT  DAMAGE  TOLERANCE  OF 
COMPOSITE  CYLINDRICAL  PANELS 


1  .  INTRODUCTION 


BACKGROUND 

Compression  of  composite  cylindrical  panels  has  been 
studied  at  the  Air  Force  Institute  of  Technology  (AFIT) 
since  the  late  1970’ s.  Wilder  [1]  presents  a  good 
background  on  the  study  of  shell  collapse.  Becker  [2] 
investigated  the  instability  of  composite  cylindrical  panels 
at  AFIT  in  1979.  Janisse  [3]  studied  the  effects  of  surface 
imperfections  and  small  cutouts  in  1982. 

Fiber  breakage,  matrix  cracking,  and/or  delamination  of 
composite  materials  may  result  from  low  energy  impacts  even 
though  no  damage  can  be  seen  on  the  surface.  Little 
information  is  available  on  impact  damage  and  impact  damage 
tolerance  of  composite  cylindrical  panels.  Much  has  been 
written  about  impact  damage  of  flat  plates  and  their 
compressive  residual  strength  after  impact  [4-10].  In  most 
cases,  the  impacted  plates  were  cut  into  tensile  or 
compression  coupons  with  the  majority  of  the  damage 
contained  in  one  coupon.  This  was  done  since  compression 
fixtures  for  flat  plates  were  not  available.  While  testing 
this  way  gives  some  insight  into  how  impact  damage  can 
weaken  a  structure,  it  confines  all  the  damage  to  one 
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concentrated  area.  The  true  structural  reaction  to  the 
damage  is  not  shown. 

Meyer  [11]  did  investigate  low  energy  impact  of 
filament-wound  Kevlar/Epoxy  panels.  He  was  able  to 
characterize  the  damage  by  calculating  the  energy  absorbed 
by  the  panel  and  by  inspecting  cross-sections  of  the  damaged 
curved  panels.  Wilder  [1]  performed  a  study  where  impact 
damage  to  a  cylindrical  Graphite/Epoxy  panel  was  simulated 
by  implanting  a  teflon  insert  into  the  panel  when  it  was 
manufactured.  The  insert  left  a  completely  debonded  region 
between  two  plies  similar  to  the  delamination  that  occurs  in 
flat  plates. 

PURPOSE 

The  purpose  of  this  thesis  was  to  study  the  damage 
tolerance  of  composite  cylindrical  panels  subjected  to  a 
single  low  energy  impact.  Two  materials  and  three  ply 
orientations  were  studied.  The  compressive  residual 
strength  was  obtained.  Another  purpose  was  to  investigate 
and  characterize  impact  damage  of  composite  cylindrical 
panels.  Ultrasonic  C-scan,  stereo  X-ray,  and 
microphotography  were  used  to  perform  the  investigation. 

SCOPE 

One  set  of  panels  tested  in  this  thesis  was  made  of 
AS4/3501-6  Graphite/Epoxy  with  ply  orientations  of  [0/90]2s, 
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[±45/90/0]s,  and  [0/90]3s> 


Panels  with  the  same  ply 


orientations  but  made  of  Kevlar/Polyester,  which  is  a  weave 
material,  were  also  tested.  Panels  from  each  group  were 
tested  in  compression  to  establish  the  original  strength  of 
the  panels.  Similarly,  panels  from  each  group  were 
subjected  to  low  energy  impacts  of  between  1-3  ft-lbs. 
C-scans  and  stereo  X-rays  were  taken  of  these  panels  to 
determine  damage  size.  They  were  then  tested  in  compression 
to  determine  residual  strength.  Lastly,  panels  from  each 
group  were  subjected  to  higher  impact  energies  (8-17 
ft-lbs).  C-scans  and  microphotographs  were  taken  of  these 
panels  in  an  effort  to  characterize  the  damage. 
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2 .  Experimentation 

The  main  purpose  of  this  thesis  is  to  investigate  how 
well  composite  cylindrical  panels  made  of  Graphite/Epoxy  and 
Kevlar/Polyester  resist  low  velocity  impacts.  Testing  this 
resistance  took  the  form  of  finding  the  residual  strength  of 
the  panels  after  they  have  been  impacted.  With  this  in 
mind,  panels  of  each  material  were  subjected  to  low  velocity 
impacts  and  tested  in  compression.  In  addition  to  residual 
strength  testing,  several  techniques  were  used  to  capture 
and  characterize  the  damage  from  a  material  point  of  view. 

Impact  of  Composite  Cvl indrical  Panels 

A  great  deal  of  work  has  been  done  in  the  area  of 
impact  on  composite  flat  plates  [4-10].  However,  the  author 
could  find  very  little  published  work  regarding  impact  of 
composite  cylindrical  panels  (or  shells  for  that  matter). 
While  it  is  not  the  primary  purpose  of  this  thesis  to  study 
the  dynamics  of  impact,  a  brief  discussion  on  the  phenomenon 
as  it  relates  to  the  testing  of  the  cylindrical  panels  is 
presented.  The  approach  is  from  the  impact  testing  point  of 
view. 

Impact  tests  were  carried  out  on  Graphite/Epoxy  and 
Kevlar/Pol yester  panels  with  the  geometry  shown  in  Figure 
2.1.  Displacements  in  the  x-  and  z-directions  are  measured 
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Figure  2.1  -  Panel  Geometry 
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in  units  of  length  while  displacements  in  the  y-direction 
are  measured  in  radians. 

The  panels  were  subjected  to  a  single,  transverse  (-z 
direction)  impact  normal  to  the  surface  of  the  panels. 

A  GRC  8250  Dynatup  Impact  Test  Machine  (hereafter  called 
Oynatup)  made  by  General  Research  Corporation,  was  the  test 
fixture  used  (see  Figure  2.2)  [12].  The  impactor  used  was  a 
1/2"  diameter  tup  weighing  approximately  5.5  pounds  and 
instrumented  with  an  accelerometer  to  measure  the  contact 
force  between  the  tup  and  the  panel.  The  data  obtained  from 
the  accelerometer  is  acquired  and  analyzed  by  a  GRC  730-1 
Instrumented  Impact  Test  Data  System. 

Two  quantities  important  in  discussing  impact  problems 
are  the  impact  energy  and  the  impact  velocity.  The  Dynatup 
Test  Data  System  calculates  the  impact  energy  from  the 
simple  relation 


E,  =  Wh  (2.1  ) 

where  W  is  the  weight  of  the  impactor  in  pounds  and  h  is  the 
drop  height  in  feet.  The  impact  energy  is  also  simply  the 
kinetic  energy  of  the  impactor  at  the  instant  it  contacts 
the  panel . 


Ej  =  1/2(mV2) 


(2.2) 
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Realizing  that  the  impactor  weight  is  just  W=mg,  we  can 
substitute  this  into  Equation  (2.1),  set  Equations  (2.1)  and 
(2.2)  equal  to  each  other  and  solve  for  the  impact  velocity, 
V. 


V 


(2.3) 


where  g  is  the  local  acceleration  of  gravity  (nominally 
32.174  ft/sec2). 

The  Dynatup  actually  calculates  the  impact  velocity  by 
use  of  a  velocity  flag.  A  metal  ‘‘flag"  similar  to  the  one 
shown  in  Figure  2.3  is  attached  to  the  impactor.  The  flag 
is  adjusted  so  that  just  prior  to  impact  it  trips  a  beam  of 
light.  The  beam  remains  off  as  the  first  prong  passes 
through  the  beam.  After  the  first  prong  passes  through  the 
beam,  the  beam  is  then  intact  as  the  space  between  the 
prongs  pass  through.  The  second  prong  then  breaks  the  beam 
again.  The  data  acquisition  system  measures  the  amount  of 
time  that  passes  between  the  two  breaks  in  the  light  beam. 
Knowing  the  distance  between  the  edges  of  the  two  prongs, 
the  system  calculates  the  impact  velocity  by  dividing  this 
distance  by  the  measured  time. 

Information  available  from  the  Dynatup  on  an  impact 
test  include  Load  vs.  Time,  Deflection  vs.  Time,  Velocity 
vs.  Time,  Load  vs.  Deflection,  Energy  vs.  Time,  or  any  other 
combination  of  these  parameters.  An  example  of  Load  vs. 
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Time  and  Energy  vs.  Time  is  shown  in  Figure  2.4.  The  load 
is  nothing  more  than  the  acceleration  of  the  tup  as  obtained 
from  the  accelerometer ,  multiplied  by  the  the  mass  of  the 
impactor.  The  Deflection  vs.  Time  curve  shown  in  Figure  2.5 
is  obtained  by  integrating  the  acceleration  vs.  time  curve 
twice  with  respect  to  time.  (Recall  from  simple  dynamics 
that  the  acceleration  is  the  second  derivative  of 
displacement  with  respect  to  time).  Knowing  the  load  and 
deflection  at  every  instant  of  time,  a  Load  vs.  Deflection 
curve,  as  shown  in  Figure  2.6,  can  be  constructed.  Lastly, 
the  Energy  vs.  Time  curve  shown  in  Figure  2.4  is  found  by 
integrating  the  Load  vs.  Deflection  curve  [13]  since 

E  =  J  F  dx  (2.4) 

One  last  thing  to  discuss  is  the  kinematic  relations 
for  a  cylindrical  panel  and  how  they  might  affect  impact 
test  (and  damage)  results  as  compared  to  flat  plates.  Even 
though  the  structure  under  study  here  is  a  cylindrical 
panel,  the  kinematics  are  the  same  as  for  a  complete 
cylindrical  shell.  Thus,  the  nonlinear  strains  and 
curvatures  for  a  cylindrical  shell  as  a  function 
displacements  can  be  written  as  [14] 
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Figure  2.4 

Example  Plot  of  Load  and  Energy  vs.  Time  for  Impact  Test 
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Example  Plot  of  Deflection  vs.  Time  for  Impact  Test 
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Example  Plot  of  Load  vs.  Deflection  for  Impact  Test 


2-10 


where 


e  = 


e  = 


u,  +  1/2  p~ 
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xy 
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x  xy 
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y  y ,  y 

K  =1/2f0  /  R  +  P  } 

xy  ^  >‘>y  y,=<  ; 


(2.5) 


(2.6) 


R  Is  the  radius  of  curvature  of  the  cylindrical  shell  and 
u,  v,  and  w  are  the  middle  surface  displacements.  Equations 
(2.5)  and  (2.6)  represent  the  Sanders  kinematic  relations 
for  shells  specialized  for  cylindrical  coordinates  [14]. 
These  equations  apply  under  the  following  assumptions  [14]: 

a.  The  shell  is  thin,  that  is  the  thickness-to-radius 
ratio  is  much  less  than  1  (h/R  <<  1). 

b.  The  strains  are  small  compared  to  unity. 

c.  Transverse  shear  and  normal  stresses  are  negligible 
(plane  stress  condition). 
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d.  The  middle  surface  strains  are  of  the  order  of 
magnitude  of  the  total  strain  squared. 

e.  The  rotations  0  and  I?  are  of  the  order  of 

x  y 

magnitude  of  the  total  strain. 

f.  The  rotations  relative  to  the  normal  to  the  middle 
surface  are  relatively  small. 

If  one  examines  Equations  (2.5)  and  (2.6),  one  sees 
that  there  are  additional  terms  in  the  strain  and  curvature 
equations  for  a  cylindrical  shell  brought  about  by  the 
curvature  of  the  shell.  Thus,  the  curvature  provides  some 
additional  strain  energy  to  the  structure.  Given  this,  one 
would  expect  that  a  greater  amount  of  force  (or  energy) 
would  be  required  to  deflect  a  curved  shell  the  same  amount 
as  a  flat  plate.  Or,  looking  at  it  from  an  impact  damage 
point  of  view,  it  would  require  more  force  (or  energy)  to 
damage  a  cylindrical  shell  than  it  would  to  damage  a  flat 
plate  of  the  same  thickness  and  ply  orientation.  The  extra 
energy  is  required  to  overcome  the  additional  strain  energy 
inherent  in  the  structure. 

With  this  brief  introduction  to  impact  and  how  the 
Dynatup  calculates  its  results,  we  can  proceed  with  a 
discussion  of  the  impact  testing. 

Impact  Testing  Procedure 

Since  the  Dynatup  could  not  accommodate  curved  panels, 
curved  aluminum  blocks  had  to  be  manufactured  to  support  the 


2-12 


panels  during  the  impact  tests.  Figure  2.7  shows  one  of 
these  blocks  with  a  Graphite/Epoxy  panel  resting  on  it.  A 
steel  cover  plate,  also  shown  in  Figure  2.7,  was  made  to 
hold  down  the  panels  during  the  tests.  A  5"  x  5”  cutout  was 
made  in  both  the  cover  plate  and  the  aluminum  support  blocks 
to  allow  the  impactor  to  strike  the  panel  and  to  allow  the 
panel  to  deflect.  The  clamped  boundary  conditions  are 
simulated  at  all  four  edges.  The  panel  is  centered  on  the 
block  such  that  the  center  of  the  panel  is  aligned  with  the 
centers  of  the  cutouts.  The  cover  plate  is  then  bolted  into 
place  with  its  cutout  aligned  with  the  cutout  of  the  support 
block.  The  final  assembly  of  the  panel  on  its  support  block 
is  shown  in  Figure  2.8. 

The  panel -support  block  assembly  is  placed  on  an 
adjustable  support  plate  underneath  the  impactor.  As  was 
mentioned  earlier,  the  mass  of  the  impactor  used  was 
approximately  5.5  lbs.  Therefore,  the  impact  velocity  and 
subsequently  the  energy  is  varied  by  changing  the  drop 
height.  For  the  Dynatup,  the  drop  height  was  set  using  a 
tape  measure.  The  height  can  be  adjusted  by  using  an 
Up-Down  switch  on  the  control  panel  shown  in  Figure  2.2. 

A  more  accurate  way  to  set  the  drop  height  is  to 
calculate  the  impact  velocity  associated  with  the  desired 
drop  height  using  Equation  (2.3).  Then,  placing  a  scrap 
panel  on  the  support  block,  adjust  the  height  of  the 
impactor  until  the  velocity  determined  by  performing  a 
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Figure  2.8  -  Assembled  Cover  P 


velocity  check  [13]  matches  the  calculated  velocity.  This 
method,  however,  is  still  not  very  easy  to  do  on  the  Dynatup 
because  there  is  no  fine-tune  or  vernier  adjustment 
capability  on  the  drop  height. 

Once  the  drop  height  is  set,  a  brake  is  set  to  catch 
the  impactor  as  it  rebounds  off  the  panel,  preventing 
multiple  impacts.  The  impact  test  is  ready  to  begin.  The 
data  acquisition  system  is  started  [13],  a  switch  is  thrown 
to  arm  the  system,  and  a  release  button  is  pressed  to 
electronically  release  the  impactor.  The  test  data  is 
acquired  and  stored  to  disk  so  that  at  any  time  after  the 
test,  the  data  may  be  analyzed. 

Minimum  Damage  Energy 

O’ Kane  and  Benham  [5]  showed  in  their  work  that  there 
is  an  impact  energy  level  below  which  no  damage  is  sustained 
by  a  composite  flat  plate.  Since  the  goal  is  to  test  panels 
for  residual  strength,  the  panels  must  be  impacted  with 
enough  energy  (or  force)  to  cause  damage.  With  this  in 
mind,  the  minimum  energy  levels  required  to  cause  damage  to 
Graphite/Epoxy  panels  with  ply  orientations  of  [0/90]^, 
[±45/90/0] s ,  and  [0/90]gs  were  determined. 

A  trial  and  error  method  was  used  to  determine  the 
minimum  energy  levels.  The  [0/90]2s  and  [±45/90/0]s 
orientations  were  initially  impacted  with  a  velocity  of  9.67 
ft/sec.  The  [0/90]3g  was  initially  impacted  with  a  velocity 
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of  11.85  ft/sec.  These  initial  impacts  were  subjected  to 
ultrasonic  C-scan  to  determine  the  extent  of  the  damage,  if 
any.  Based  on  the  results  of  the  C-scans,  the  impact 
velocity  (energy)  was  increased  or  decreased  accordingly. 
This  process  was  repeated  until  the  minimum  energy  level 
that  caused  damage  was  determined.  Table  2.1  summarizes  the 
minimum  damage  energies  for  the  three  types  of 
Graphite/Epoxy  panels. 


Table  2.1 

Minimum  Damage  Energy  for  Gr/Ep  Panels 


C0/90] gs 

[±45/90/0]s 

[0/90] 3s 

Energy  (ft- lbs) 

1.27 

0.99 

2.46 

By  looking  at  Table  2.1,  we  see  how  the  thickness  and 
ply  orientations  affect  the  panels  ability  to  absorb  the 
impact.  The  [±45/90/0]  panel  was  the  most  easily  damaged, 

S 

while  the  thicker  [0/90]3s  panel  required  more  energy  to 
damage  it.  A  more  in  depth  discussion  of  the  impact  damage 
will  be  presented  in  Chapter  3. 

In  order  to  insure  that  enough  damage  is  present  in  the 
panel  to  see  a  difference  in  the  residual  strength  tests, 
the  impact  energies  used  in  the  actual  testing  were 
increased  20%  above  those  shown  in  Table  2.1.  These 
increased  energy  levels  are  shown  in  Table  2.2. 


Table  2.2 

Energy  Levels  Used  for  Impact  Tests 


[0/90]2g 

[±45/90/0]g 

[0/90]3s 

Energy  (ft-lbs) 

1 .52 

1.19 

2.92 

Impact  Tests 

Using  the  impact  energies  shown  in  Table  2.2  and 
following  the  impact  testing  procedures  described  earlier, 
three  Graphite/Epoxy  and  three  Kevl ar/Pol yester  panels  were 
impacted  with  the  same  energy  for  the  given  ply  orientation. 
For  example,  three  [0/90]2g  panels  of  each  material  were 
impacted  with  1.52  ft-lbs  of  energy. 

Figures  2  ^  '.14  show  one  Load  and  Energy  vs.  Time 
plot  for  eacn  ^anel  type.  The  small  oscillations  apparent 
in  the  Load  vs.  Time  curve  result  from  noise  in  the  Dynatup 
as  thb  impactor  falls,  strikes  the  panel,  and  rebounds.  In 
the  Graphite/Epoxy  panels  (Figures  2.9,  2.11,  and  2.13),  we 
notice  larger  oscillations  in  the  force,  corresponding  to 
the  panels  being  damaged.  The  Kevlar/Polyester  panels  hit 
with  the  same  velocities  show  no  such  indications  of  damage. 
However,  C-Scans  indicate  some  damage  resulted  from  the 
impact. 

For  each  ply  orientation,  we  also  notice  the  impact 
generates  a  greater  force  in  less  time  in  the  Graphite/Epoxy 
panels  as  compared  to  the  Kevlar/Polyester  panels.  For 
instance,  the  maximum  load  for  [0/90]2g  Graphite/Epoxy  is 
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Figure  2.10 

Load  and  Energy  vs.  Time  for  [0/90],  Kev/Poly  Panels 
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Figure  2.12 

Load  and  Energy  vs.  Time  for  [±45/90/01  Kev/Poly  Panels 
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Figure  2.13 

Load  and  Energy  vs.  Time  for  C0/90].  Gr/Ep  Panels 
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199.3  lbs  and  occurs  at  7.02  msec.  For  a  [0/90] _ 
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Kevlar/Pol yester  panel,  the  maximum  load  is  157.8  lbs 
occurring  at  9.13  msec.  The  Kevlar/Polyester  panel 
dissipates  some  of  the  energy  by  deflecting  more  because  it 
is  more  flexible.  This  greater  deflection  shows  up  in  the 
greater  amount  of  time  required  for  the  panel  to  reach  the 
peak  load. 

Compression  Testing  Procedure 

The  test  fixture  used  to  compress  the  panels,  shown  in 
Figure  2.15,  was  provided  by  the  Air  Force  Flight  Dynamics 
Laboratory.  The  panel  is  placed  into  the  fixture  between 
the  top  and  bottom  plates  with  the  top  and  bottom  curved 
edges  of  the  panel  clamped.  The  vertical  edges  are 
restrained  by  a  simple  support  knife  edge.  The  clamped  edge 
condition  is  achieved  by  placing  a  series  of  six  1/2"  steel 
chucks  along  the  top  and  bottom  edges.  The  edges  of  the 
chucks  have  a  12"  radius  of  curvature  so  that  they  butt 
evenly  against  the  panel.  Each  chuck  is  held  in  place  by 
two  set  screws.  See  Reference  [1]  for  a  more  complete 
description  of  the  compression  fixture.  Mathematically,  the 
boundary  conditions  are 

Top  edge:  u=free;  v=w=w, x=w, y=0 

Bottom  edge:  u=v=w=w,  =w,  =0 

x  y 

Vertical  edges:  u=w ,  =free;  v=w=w,  =0 

x  y 
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Figure  2.16  -  Strain  Gauge  Placement 
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As  the  panel  is  compressed,  a  load  cell  measures  the 
total  load  applied  to  the  top  edge  while  a  LVDT  (Linear 
Variable  Differential  Transformer)  measures  the  vertical 
axial  displacement,  u.  Additionally,  unidirectional  strain 
gauges  were  mounted  on  each  panel  as  shown  in  Figure  2.16. 
The  two  strain  gauges  near  the  top  edge  of  the  panel  provide 
information  to  insure  the  load  is  being  applied  as  evenly  as 
possible  across  the  top  edge.  The  two  strain  gauges  in  the 
center  of  the  panel  (placed  back  to  back)  measure  the  strain 
reversal  that  occurs  when  the  panel  buckles. 

After  a  panel  is  placed  in  the  fixture,  the  1/2" 
chucks  and  vertical  knife  edges  are  loosely  set  in  place.  A 
preload  of  approximately  200  lbs  for  Graphite/Epoxy  panels 
(100  lbs  for  Kevlar/Polyester  panels)  is  applied.  The 
preload  makes  sure  the  top  and  bottom  edges  of  the  panel  are 
flush  against  the  base  plates.  The  support  bolts  and  set 
screws  are  then  tightened  finger  tight.  The  preload  is  then 
eased  to  approximately  100  lbs  for  Graphite/Epoxy  panels  (50 
lbs  for  Kevlar/Polyester  panels)  and  the  support  bolts  and 
set  screws  tightened  with  a  wrench  and  screwdriver.  This 
preload  procedure  prevents  movement  of  the  panel  before  it 
begins  to  take  on  the  load. 

Once  the  preload  procedure  is  complete,  the  panels  are 
loaded  in  compression.  A  uniform  displacement  of  0.05 
inches/minute  is  applied  to  the  top  edge  with  a  load  cell 
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measuring  the  resulting  force  and  the  vertical  LVDT 
measuring  the  top  edge  displacement. 

The  load  is  applied  until  the  panels  buckle.  For  the 
Graphite/Epoxy  panels,  a  loud  "pop"  can  be  heard  when  the 
panels  snap  into  their  buckled  shape.  For  the 
Kevlar/Polyester  panels,  no  “pop"  was  heard  due  to  the 
greater  flexibility  of  the  material.  The  movement  of  the 
panel  as  it  buckles  is  not  as  quick,  thus  it  becomes  more 
difficult  to  distinguish  the  buckling  load.  As  a  result, 
the  output  from  the  load  cell  was  monitored  and  the  test 
stopped  when  the  load  began  to  decrease.  The  maximum  load 
observed  was  called  the  buckling  load.  Once  the  buckling 
load  was  determined,  the  test  was  stopped,  the  load  removed 
and  the  panel  removed  from  the  test  fixture. 

Compression  of  Panels 

Three  undamaged  Graphite/Epoxy  and  three  undamaged 
Kevlar/Polyester  panels  of  each  ply  orientation  were 
compressed  to  establish  the  baseline  strength  of  the  panels 
The  loads  measured  will  be  used  for  comparison  purposes  to 
determine  how  the  impact  damage  weakens  the  panels.  Table 
2.3  summarizes  the  collapse  loads  of  the  undamaged  panels. 
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Table  2.3 

Collapse  Loads  of  Undamaged  Panels 


[0/90]2s 

[±45/90/0]s 

ro/9°]3s 

Load 

(lbs) 

Gr/Ep 

2802 

3551 

7667 

Kev/Pol 

813 

1103 

2408 

Three  Graphite/Epoxy  and  three  Kevlar/Polyester  panels 
of  each  ply  orientation  were  impacted  with  the  energies 
shown  in  Table  2.2  using  the  procedures  described  earlier. 
These  panels  were  then  compressed  to  determine  the  residual 
strength  of  the  panels.  The  residual  strengths  are 
summarized  in  Table  2.4. 


Table  2.4 

Collapse  Load  of  Damaged  Panels 


C0/90]2s 

[±45/90/0] 

t 0/9°) 3s 

Load 

(lbs) 

Gr/Ep 

2750 

3359 

7492 

Kev/Pol 

769 

1 118 

2346 

Comparing  the  values  in  Tables  2.3  and  2.4,  the  damaged 
panels  were  2-6%  weaker  than  the  undamaged  panels.  The  only 
exception  is  the  [±45/90/0]  Kevlar/Polyester  panels.  The 

S 

average  strength  of  the  damaged  panels  was  15  lbs  greater 
than  that  of  the  undamaged  panels.  One  should  not  infer 
from  this  that  this  particular  panel  becomes  stronger  after 
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it  is  impacted.  Recall,  the  values  in  Tables  2.3  and  2.4 
were  based  on  a  sample  of  three  panels.  One  must  also  keep 
in  mind  the  relative  insensitivity  of  the  compression 
fixture  alluded  to  earlier.  And,  there  is  always  a  certain 
amount  of  experimental  scatter  to  live  with.  Thus,  the 
results  of  this  test  for  [±45/90/0]  Kevlar/Pol yester  panels 

S 

should  be  interpreted  as  the  impact  energy  used  having 
little  or  no  effect  on  the  panel’s  compressive  strength. 
However,  ultrasonic  C-scans  of  the  damaged  panels  taken 
after  the  one  time  compression  load  showed  the  damaged  area 
in  the  panel  had  spread.  The  damaged  area  before  the 
compression  test  was  approximately  1/2"  in  diameter.  After 
the  compression  testing  of  the  damaged  panels,  the  damaged 
area  of  the  Graphite/Epoxy  panels  did  not  grow  appreciably, 
while  the  damaged  area  of  the  Kevlar/Pol yester  panels  grew 
to  approximately  1-1/4"  to  2"  in  diameter,  depending  on  the 
ply  orientation.  It  would  be  reasonable  to  assume  that  if 
the  panels  were  loaded  in  compression  again,  they  would 
buckle  at  lower  loads  than  shown  in  Table  2.4. 

The  results  of  this  experimentation  tend  to  verify 
Wilder’s  technique  [1]  of  implanting  teflon  inserts  to 
simulate  impact  damage.  His  results  showed  a  decrease  in 
global  strength  of  13. 2%  for  a  quasi-isotropic  ply 
orientation.  One  must  bear  in  mind  that  Wilder’s  technique 
would  only  give  one  an  idea  of  the  global  behavior  of  the 
panel.  Near  the  impacted  region,  the  behavior  may  be  quite 
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different.  In  an  impacted  panel,  the  damaged  area  may  still 
have  some  strength  left  in  it  depending  on  the  amount  of 
energy  imparted  to  the  panel.  The  teflon  inserts  Wilder 
used  provides  a  totally  debonded  region.  Thus,  the  teflon 
insert  technique  would  tend  to  yield  a  conservative  estimate 
of  an  impact  damaged  panel’s  global  response  to  compression. 

Figures  2.17-2.22  show  plots  comparing  load  vs. 
in-plane  displacement,  u,  for  damaged  and  undamaged  panels 
of  each  material  and  ply  orientation  used  to  generate  Tables 
2.3  and  2.4.  In  each  case,  little  difference  in  the 
stiffness  of  the  damaged  panels  can  be  seen  compared  to  the 
undamaged  panels.  Only  the  collapse  load  differs. 
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Load  (lbs) 


Load  (lbs) 


Displacement,  U  (*10**-3  in.) 


Figure  2.19 

Load  vs.  Displacement  for  [±45/90/0]  Gr/Ep 

s 

Damaged  and  Undamaged  Panels 
(Impact  Energy  =  1.19  ft- lb) 
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Load  (lbs) 


Load  (lbs) 


Figure  2.22 

Load  vs.  Displacement  for  [0/90] 3s  Kev/Poly 

Damaged  and  Undamaged  Panels 
(Impact  Energy  =  2.92  ft-lb) 
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3.  DAMAGE  CHARACTERIZATION 

In  this  chapter,  we  take  a  closer  look  at  the  damage 

sustained  by  the  panels  from  a  low  velocity  impact.  First, 

stereo  X-rays  were  taken  of  the  Graphite/Epoxy  panels  used 

to  determine  the  minimum  impact  energy  to  cause  damage. 

Second,  a  series  of  photographs  and  ultrasonic  C-Scans  were 

taken  of  Graphite/Epoxy  and  Kevlar/Polyester  panels  impacted 

with  velocities  ranging  from  9.62-13.89  ft/sec  (energies  of 

7.96-16.57  ft-lb).  We  therefore  get  an  idea  of  what  happens 

as  the  impact  energy  increases.  Lastly,  [0/90]2g 

Graphite/Epoxy,  [0/90]2g  Kevlar/Pol yester ,  [±45/90/0]g 

Graphite/Epoxy,  and  [±45/90/0]_  Kevlar/Pol yester  panels  were 

s 

cross-sectioned  through  the  impact  point  and  viewed  under  a 
microscope.  This  shows  the  damage  through  the  thickness  of 
the  panel . 

Stereo  X-rav  of  Panels 

Stereo  X-rays  provide  a  way  of  getting 
three-dimensional  information  on  a  damaged  composite  by 
combining  two  images  taken  at  slightly  different  angles  into 
one  apparent  three-dimensional  image.  The  three-dimensional 
image  is  produced  much  the  same  way  the  human  eye  perceives 
it.  Each  eye  sees  an  image  at  a  slightly  different  angle 
due  to  the  separation  between  the  eyes.  The  brain  then 
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merges  the  two  images  together  to  form  a  three-dimensional 
image. 

For  this  thesis,  the  stereo  X-rays  were  prepared  as 
follows  [15].  A  1/8"  hole  was  drilled  in  the  damaged  area 
of  the  panel  to  allow  a  dye  penetrant,  tetrabromoethane 
(TBE),  to  soak  into  the  damaged  area  via  capillary  action. 
The  TBE  provides  a  color  contrast  between  the  damaged  and 
undamaged  areas  of  the  panel.  Next,  film  is  placed  behind 
the  panel  and  exposed  by  using  an  X-ray  tube.  Two  exposures 
are  taken  at  ±15°on  either  side  of  a  plane  perpendicular  to 
the  panel  (see  Figure  3.1).  These  two  images  may  be 
recombined  by  using  a  "stereo  viewer"  to  obtain  a 
three-dimensional  view  of  the  damaged  area. 

Figures  3. 2-3. 4  are  photographs  of  the  X-rays  for  the 
minimum  damage  of  each  type  of  Graph ite/Epoxy  panel.  In 
each  figure,  some  damage  to  the  fibers  as  well  as  some 
cracking  of  the  fibers  and  matrix  can  be  seen.  The  figures 
also  show  the  extent  and  shape  of  the  delamination  that 
occurs.  It  is  not  very  easy  to  pick  out  of  the  photographs 
by  eye,  but  under  the  stereo  viewer,  one  can  see  the 
delamination  occurring  at  several  different  layers  through 
the  thickness  of  the  panel.  We  only  get  a  feel  of  relative 
positioning  though.  It  is  still  difficult  to  pinpoint 
between  which  two  plies  the  delaminations  actually  occur. 
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Figure  3.1  -  Stereo  X-Ray  Geometry 
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Figure  3.2 

Stereo  X-Ray  of  [0/90 ]2  Gr/Ep  Impacted  Panel 
(Impact  Energy  =  1.27  ft-lb) 
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Figure  3.3 

Stereo  X-Ray  of  [±45/90/01  Gr/Ep  Impacted  Panel 

© 

(Impact  Energy  =  0.99  ft- lb) 
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Figure  3.4 

Stereo  X-Ray  of  C0/90]3s  Gr/Ep  Impacted  Panel 
(Impact  Energy  =  2.46  ft-lb) 


Photographs  and  C-Scans  of  Panels 


In  this  phase  of  testing,  Graphite/Epoxy  and 
Kevlar/Polyester  panels  were  impacted  with  increasing 
velocities  and  C-Scanned.  By  doing  this,  we  can  track  how 
the  damage  spreads  as  the  energy  imparted  to  the  panel 
increases. 

We  begin  by  considering  [0/90]2s  Graphite/Epoxy  panels 
impacted  with  velocities  ranging  from  9.62-12.68  ft/sec. 
Figures  3.5-3. 7  shows  the  impacted  side,  a  C-Scan  of  the 
internal  damage,  and  the  back  side  of  Graphite/Epoxy  panels 
impacted  at  9.62  ft/sec,  11.34  ft/sec,  and  12.68  ft/sec 
respectively.  As  seen  from  the  C-Scans,  the  internal  damage 
increases  as  the  impact  velocity  increases.  The  damage 
occurs  in  the  form  of  a  circular  region  at  the  point  of 
impact  and  is  approximately  the  same  diameter  as  the 
impactor  (1/2").  Additionally,  two  triangular  regions  are 
damaged  on  each  side  of  the  impact  point  in  the 
circumferential  direction.  Since  the  impactor  is 
hemispherical,  one  would  not  expect  any  directional 
preference  of  the  damage.  However,  the  curvature  of  the 
panel  seems  to  dictate  the  direction  of  the  damage.  This 
also  helps  to  explain  the  fiber  splitting  and  fracture  that 
occurs  on  the  back  side  of  the  panel.  Locally,  the  panel 
undergoes  a  high  concentration  of  bending  and  tension.  The 
panel  is  able  to  deflect  ci rcumferential ly ,  but  since  it  is 
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Figure  3.5  -  Impacted  [0/90] 2g  Gr/Ep  Panel 

(Velocity  =  9.62  ft/sec;  Impact  Energy  =  7.96  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


Figure  3.6  -  Impacted  C0/90]28  Gr/Ep  Panel 

(Velocity  =  11.34  ft/sec;  Impact  Energy  =  11.06  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


3-9 


Figure  3.7  -  Impacted  f0/90]2s  Gr/Ep  Panel 

(Velocity  =  12.08  ft/ sec;  Impact  Energy  =  13.77  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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much  stiffer  longitudinally,  the  fibers  break  in  response  to 
the  contact  force  created  by  the  impact. 

In  Figures  3.8-3.10  we  look  at  [0/90]2s 
Kevlar/Polyester  panels  impacted  with  the  same  velocities  as 
the  equivalent  Graphite/Epoxy  panels.  The  C-Scans  show  the 
damage  growing  in  the  ci rcumferential  direction  as  the 
impact  velocity  increases,  just  as  for  Graphite/Epoxy.  The 
width  of  the  damaged  area  remains  approximately  the  diameter 
of  the  impactor.  We  do  notice  something  different  about  the 
Kevlar/Polyester  panels  however.  As  the  velocity  increases 
from  11.34  ft/sec  to  12.68  ft/sec,  the  shape  of  the  internal 
damage  becomes  more  circular.  It  appears  as  if  the  type  or 
shape  of  damage  that  occurs  depends  on  the  impact  velocity 
(or  possibly  energy).  As  the  velocity  increases,  the  damage 
to  the  back  side  changes  from  a  slit  where  the  fibers 
fracture  to  a  cross  where  the  fibers  fracture  in  both 
directions.  Also,  we  notice  some  damage  occurring  along 
diagonals  extending  from  the  impact  point  to  the  corners 
where  the  panel  is  supported  by  the  support  block  and  cover 
plate.  A  permanent  indentation  remains  as  a  result  of  the 
impact.  It  appears  that  as  a  result  of  the  transverse 
displacement,  the  material  stretched,  causing  damage  along 
the  lines  described. 

Figure  3.11  shows  a  [±45/90/0]  Graphite/Epoxy  panel 

S 

impacted  with  a  velocity  of  9.82  ft/sec.  Only  one  panel  of 
this  type  was  available  for  testing,  therefore  no  real 
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Figure  3.8  -  Impacted  [0/90]2s  Kev/Poly  Panel 

(Velocity  =  9.62  ft/sec;  Impact  Energy  =  7.96  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


Figure  3.9  -  Impacted  [0/90]2s  Kev/Poly  Panel 

(Velocity  =  11.34  ft/sec;  Impact  Energy  =  11.06  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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conclusions  can  be  drawn  by  studying  it  alone.  However,  we 
can  compare  it  to  a  Kevlar/Polyester  panel  to  see  how  the 
two  materials  respond  for  the  same  ply  orientation  and 
impact  velocity.  The  impacted  side  does  not  appear  to  be 
too  badly  damaged,  however,  on  the  back  side,  we  see  a  great 
deal  of  fiber  splitting  and  breakage.  From  the  C-Scan,  most 
of  the  internal  damage  appears  to  spread  along  the  -45° 
direction.  Comparing  this  panel  to  the  equivalent 
Kevlar/Polyester  panel  shown  in  Figure  3.12,  we  see  similar 
internal  damage  from  the  C-Scan.  The  damage  is  more 
circular  and  seems  to  show  a  tendency  to  spread  in  the  ±45° 
directions.  Again,  the  damage  extending  from  the  corners  to 
the  impact  point  is  present. 

Lastly,  we  consider  [0/90]3s  panels.  Figures  3.13-3.16 
illustrate  the  damage  of  Graphite/Epoxy  panels  where  Figures 
3.17-3.20  show  damaged  Kevl ar/Pol yester  panels.  Both  sets 
of  panels  were  impacted  with  velocities  ranging  from 
11.75-13.89  ft/sec. 

For  the  Graphite/Epoxy  panels,  we  notice  very  little 
damage  on  the  impacted  side.  This  contrasts  with  the 
[0/90] 2S  panels  which  showed  a  little  more  damage  on  the 
impacted  surface.  The  back  sides  however  show  the  exact 
same  features  as  the  thinner  panels  (fiber  splitting  and 
breakage).  The  internal  damage  (delamination)  evidenced 
from  the  C-Scans  are  somewhat  different  than  for  the 
[0/90]2s  panels.  While  the  shapes  of  the  damaged  areas 
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Figure  3.11  -  Impacted  [±45/90/0]g  Gr/Ep  Panel 

(Velocity  =  9.82  ft/sec;  Impact  Energy  =  8.25  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


Figure  3.12  -  Impacted  [±45/90/0]s  Kev/Poly  Panel 

(Velocity  =  9.82  ft/sec;  Impact  Energy  =  8.25  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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Figure  3.13  -  Impacted  [0/90]3s  Gr/Ep  Panel 

(Velocity  =  11.75  ft/sec;  Impact  Energy  =  11.81  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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Figure  3.14  -  Impacted  [0/90]3s  Gr/Ep  Panel 

(Velocity  =  12.03  ft/sec;  Impact  Energy  =  12.41  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  8ack  Side 
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Figure  3.15  -  Impacted  f0/90]3s  Gr/Ep  Panel 

(Velocity  =  12.68  ft/sec;  Impact  Energy  =  13.77  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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Figure  3.16  -  Impacted  [0/90]^  Gr/Ep  Panel 

Velocity  =  13.89  ft/sec;  Impact  Energy  =  16.57  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


Figure  3.18  -  Impacted  [0/90]3g  Kev/Poly  Panel 

(Velocity  =  12.03  ft/sec;  Impact  Energy  =  12.41  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 


(Velocity  =  12.68  ft/sec;  Impact  Energy  =  13.77  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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Figure  3.20  -  Impacted  [0/90]3g  Kev/Poly  Panel 

(Velocity  =  13.89  ft/sec;  Impact  Energy  =  16.57  ft-lb) 
(a)  Impacted  Side  (b)  C-Scan  (c)  Back  Side 
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appear  to  be  similar,  it  is  surprising  to  note  that  the 
damage  area  did  not  increase  very  much  (if  at  all)  as  the 
impact  velocity  increased.  This  leads  one  to  suspect  that 
there  may  be  a  range  of  impact  energies  (or  velocities) 
where  the  amount  of  damage  levels  off.  Indeed,  O’Kane  and 
Benham  [5]  plotted  experimentally  obtained  damage  areas  as  a 
function  of  impact  energy.  For  several  of  the  ply 
orientations  they  studied,  the  curves  appear  to  flatten  out 
indicating  that  the  damage  area  becomes  more  or  less 
constant  over  a  range  of  impact  energies.  The  same  features 
are  noticed  for  [0/90]3s  panels  made  of  Kevl ar/Pol yester 
shown  in  Figures  3.17-3.20.  The  damage  does  not  appear  to 
increase  as  the  impact  velocity  increases. 

To  briefly  summarize,  damage  to  Graphite/Epoxy  panels 

takes  the  form  of  fiber  splitting,  fiber  breakage,  and 

delamination.  The  delaminations  seem  to  be  oriented  along 

the  fiber  directions.  For  the  cases  of  [0/90] _  and 

2s 

[0/90]3s  panels,  the  curvature  of  the  panel  seems  to  lend 
some  preference  as  to  which  directions  the  delaminations 
spread . 

The  damage  of  Kevl ar/Pol yester  panels  remains  much  more 
contained  and  localized  compared  to  Graphite/Epoxy  panels. 
For  Kevlar/Polyester,  we  see  fibers  breaking  in  a  tensile 
mode,  especially  on  the  back  side,  and  we  see  no  fiber 
splitting  as  in  Graphite/Epoxy,  probably  due  to  the  weave 
construction  of  the  material.  Also,  different  modes  of 
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damage  may  exist  for  Kevlar/Polyester,  based  on  the  observed 
changes  in  damage  shape  as  the  impact  velocity  increases. 

Lastly,  panels  made  of  both  materials  exhibited  the 
characteri sti c  of  a  more  or  less  constant  damage  area  over  a 
range  of  impact  energies. 

Microphotographs  of  Impacted  Panels 

The  last  part  of  investigating  the  impact  damage 
consisted  of  cross-sectioning  some  of  the  panels  and  viewing 
them  under  a  microscope.  Mi crophotographs  through  the 
impact  point  were  taken  of  these  panels  to  be  studied. 

The  following  four  panels  were  selected  to  be  studied. 

1.  [0/90]2s  Gr/Ep  (Velocity  =  12.68  ft/sec) 

2.  [0/90]2s  Kev/Poly  (Velocity  =  12.68  ft/sec) 

3.  [±45/90/0 ]  Gr/Ep  (Velocity  =  9.82  ft/sec) 

s 

4.  [±45/90/0]  Kev/Poly  (Velocity  =  9.82  ft/sec) 

s 

Metal lographic  specimens  from  each  of  these  panels  had 
to  be  prepared.  The  impacted  panels  were  first  cut  along  a 
line  through  the  impact  point  as  shown  in  Figure  3.21  using 
a  water  cooled  diamond  saw  blade.  The  water  prevented  any 
localized  heating  during  the  cutting.  A  3"  section 
surrounding  the  impact  point  was  removed  from  one  of  the 
half  panels  to  use  as  the  specimen. 
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Figure  3.21  -  Microphotograph  Specimen  Geometry 
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Each  specimen  was  placed  in  an  epoxy  resin  and  allowed 
to  set  for  approximately  48  hours.  Setting  the  specimen  in 
the  epoxy  made  polishing  easier  and  provided  stability  of 
the  specimen  for  viewing  under  the  microscope.  After  the 
epoxy  had  hardened,  each  specimen  was  ground  on  a  grinding 
wheel  until  the  cross-section  through  the  thickness  was 
exposed.  Next,  the  surface  was  sanded  using  successively 
finer  grits  of  sandpaper  (240,  320,  400,  600  grit).  Lastly, 
each  specimen  was  polished  with  an  alumina  polishing  cloth 
on  a  grinding  wheel  until  the  desired  details  could  be  seen 
under  the  microscope. 

Figures  3.22-3.25  are  mi crophotographs  showing  impact 
damage  through  the  thickness.  The  actual  width  of  the 
section  photographed  is  approximately  1/2",  the  diameter  of 
the  impactor.  The  magnification  used  to  view  the  specimens 
was  50X. 

Figure  3.22  shows  the  cross-section  of  a  [0/90]2s 
Graphite/Epoxy  panel  impacted  with  a  velocity  of  12.68 
ft/sec.  The  top  ply  of  0°  fibers  debonded  completely  from 
the  90°  ply  beneath  it.  Delaminations,  which  appear  as  dark 
horizontal  lines  in  the  photograph,  occur  between  every  ply. 
We  also  see  many  cracks  running  through  and  along  the  0° 
plies  as  well  as  numerous  small  cracks  in  the  90°  plies. 

The  horizontal  cracks  through  the  inner  0°  plies  suggest 
some  shear  failure,  possibly  resulting  from  a  shear  stress 
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o  .  The  vertical  cracks  suggest  tensile  flexural  stresses 

yz 

[6]. 

Figure  3.23  shows  the  cross-section  of  a  [0/90]2s 

Kevl ar/Pol yester  panel  impacted  by  the  same  velocity  (12.68 

ft/sec).  From  this  figure,  one  gets  a  good  feel  for  the 

weave  construction  of  Kevl ar/Pol yester .  Through  the 

thickness,  there  appears  to  be  much  less  damage  than  in  the 

equivalent  Graphite/Epoxy  panel.  The  fibers  appear  to 

remain  intact,  sustaining  no  visible  damage.  There  are, 

however,  a  number  of  cracks  in  the  matrix  material  between 

the  fibers.  Overall,  we  can  also  see  the  permanent 

deformation  of  the  impacted  region. 

Moving  on,  the  [±45/90/0]_  Graphite/Epoxy  panel 

s 

impacted  at  9.82  ft/sec  is  shown  in  Figure  3.24.  The 
cross-section  appears  to  be  literally  destroyed.  The  fibers 
are  broken  into  small  pieces  and  a  great  deal  of 
delamination  has  occurred  between  every  ply. 

Lastly,  Figure  3.25  shows  the  [±45/90/0] 

s 

Kevl ar/Pol yester  impacted  at  9.82  ft/sec.  Again,  we  see 
much  more  damage  (delamination  and  matrix  cracks)  than  for 
the  [0/90]2s  ply  orientation. 
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Figure  3.22  -  Cross-Section  of  Impacted  [0/90]  Gr/Ep  Panel 

2s 

(Impact  Velocity  =  12.68  ft/sec;  Impact  Energy  =  13.77  ft-lb) 
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Figure  3.25  -  Cross-Section  of  Impacted  [  +45/90/0]  8  Kev/Poly  Panel 
(Impact  Velocity  =  9.82  ft/sec;  impact  Energy  =  8.25  ft-lb' 


4.  ANALYSIS 


To  the  author’s  knowledge,  no  valid  model  exists  for 
predicting  the  extent  of  damage  and  corresponding  reduction 
in  strength  of  composite  cylindrical  panels  (or  shells). 

That  being  the  case,  the  analysis  performed  in  this  thesis 
concentrates  on  predicting  the  collapse  load  of  undamaged 
panels . 

An  analysis  first  performed  to  determine  the  collapse 
load  of  the  panels  assumed  a  perfect  geometry.  The 
resulting  collapse  loads  greatly  disagreed  with  the 
experimental  1 y  measured  loads  of  the  undamaged  panels. 
Depending  on  the  ply  orientation,  the  difference  between  the 
analytical  and  experimental  load  was  as  much  as  62X.  This 
prompted  an  investigation  into  geometrically  imperfect 
panels. 

Discrepancies  were  also  noted  between  the  displacements 
of  the  analytical  vs.  experimental  panels.  The  actual 
panels  were  not  as  stiff  as  the  analysis  indicates,  which  by 
itself  was  not  that  surprising.  However,  the  magnitude  of 
the  difference  was  puzzling.  As  a  result,  the  boundary 
conditions  were  relaxed  in  order  to  allow  more  structural 
flexibility.  While  relaxing  the  boundary  conditions  does 
have  an  effect  on  the  flexibility,  it  alone  does  not  account 
for  the  large  differences. 
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Nonl inear  Col  lapse  Analysis 

The  finite  element  program  STAGSC-1  was  used  to  perform 
the  analysis  of  the  composite  cylindrical  panels.  STAGSC-1 
is  an  energy  based  finite  element  code  for  performing 
structural  analysis  of  shells.  One  of  the  advantages  of 
STAGSC-1  is  its  ability  to  perform  a  nonlinear  static 
analysis.  Such  an  analysis  enabled  the  determination  of  the 
collapse  load  of  the  panels  studied  in  this  thesis.  In  the 
finite  element  formulation  of  STAGSC-1,  the  kinematic 
relations  of  a  cylindrical  shell  described  by  Equations 
(2.5)  and  (2.6)  are  replaced  with  the  nonlinear  kinematic 
relations  of  a  flat  plate.  This  is  done  because  STAGSC-1 
uses  flat  plate  elements  to  approximate  a  shell.  While  the 
thin  shell  and  small  strains  assumptions  hold,  the  inclusion 
of  the  higher  order  terms  allows  for  the  large  rotations 
that  can  occur  when  a  shell  deforms.  A  complete  derivation 
of  the  relations  is  given  in  [16].  The  middle  surface 
strains  used  by  STAGSC-1  are 
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(4.1) 


Applying  the  Ki rchhof f-Love  hypothesis  (thin  element 
with  no  distortion  of  the  cross-section)  and  including  the 


curvatures,  the  strains  at  any  point  in  the  flat  element  may 
be  expressed  as  [17] 
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where  e°  e°,  and  7°  are  given  by  Equations  (4.1).  These 

x  y  x  y 

strain  relations  are  the  ones  used  in  the  constitutive 
relations  for  composite  materials.  Appendix  B  provides  a 
brief  discussion  of  laminated  plate  theory.  The  result  of 
an  analysis  using  these  nonlinear  kinematic  relations  is  the 
equilibrium  path  the  structure  follows  as  it  deforms.  The 
peak  load  on  the  curve  is  called  the  collapse  load. 

As  was  stated  earlier,  STAGSC-1  uses  flat  elements  to 
model  the  curved  surfaces  of  a  shell.  The  two  elements 
contained  in  STAGSC-1  primarily  used  for  this  type  of 
analysis  are  the  QUAF  410  and  QUAF  411  elements.  The  QUAF 
410  is  a  quadri lateral  element  with  three  translational  and 
three  rotational  degrees  of  freedom  per  node,  allowing  24 
degrees  of  freedom  per  element.  The  QUAF  411  element,  in 
addition  to  the  same  translations  and  rotations  in  the  410 
element,  allows  an  additional  in-plane  rotation  at  each 
corner  and  also  includes  a  translational  degree  of  freedom 
at  the  midpoint  of  each  side  of  the  element.  This  amounts 
to  a  total  of  32  degrees  of  freedom  for  the  QUAF  411 


element.  Figure  4.1  illustrates  the  two  elements  and  their 
respective  degrees  of  freedom  [22]. 


Looking  at  Figure  4.1,  we  notice  a  degree  of  freedom 
not  usually  associated  with  plate  elements.  The  normal 
rotation  B^  for  the  elements  of  Figure  4.1  is  needed  to 
account  for  the  fact  that  two  adjacent  plate  elements  meet 
at  an  angle  to  represent  a  curved  surface  [17]. 

Sol ution  Technique 

The  analysis  of  an  elastic  stability  problem  involving 
shell  structures  requires  the  solution  of  a  large  system  of 
nonlinear  equations  using  a  finite  element  or  finite 
difference  method  [18].  Bellini  and  Chulya  [19]  discuss 
the  fact  that  when  the  nonlinear  response  curve  (load  vs. 
displacement)  of  a  structure  reaches  a  turning  point,  one  of 
two  things  can  happen.  A  limit  point  is  reached  where  the 
solution  stops,  or  a  bifurcation  point  is  reached  with 
several  possible  branches. 

The  solution  technique  STAGfC-1  uses  to  solve  the 
nonlinear  equations  is  the  Riks  method  [18,20]  combined  with 
the  modified  Newton-Raphson  method  as  an 
i terative-incremental  method  [19].  The  advantage  of  this 
technique  is  the  response  curve,  which  represents 
equilibrium  states,  can  be  solved  for  past  a  peak  load.  In 
essence,  the  solution  continues  and  calculates  a 
post-buckling  equilibrium  path.  Interested  readers  are 


referred  to  References  [18,19,20]  for  In  depth  discussions 
on  Riks  method  and  nonlinear  equation  solution  techniques 
using  finite  element  methods. 


Col  lapse  of  Undamaged  Panels 

Material  properties  of  AS4/3501-6  Graphi te/Epoxy  were 
determined  by  standard  tensile  tests  [21].  [0]Q  coupons 

were  used  to  determine  E1  and  t>12  while  [ 90 1 1 6  and  [±45]2 
coupons  were  used  to  determine  E2  and  G12  respectively.  The 
material  properties  determined  for  the  analysis  of 
Graphite/Epoxy  panels  are 

E1  =  20.461  *  106  psi 
E2  =  1.34036  *  106  psi 
i'12  =  0.3131 
G12  =  0.8638  *  106  psi 

Determining  the  material  properties  of  Kevlar/Polyester 
is  not  as  straight  forward  as  for  Graphite/Epoxy  because 
Kevlar/Polyester  can  not  be  made  into  unidirectional 
coupons.  Also,  STAGSC-1  can  not  analyze  a  weave  material, 
so  an  equivalent  orthotropic  laminate  had  to  be  developed 
to  model  Kevlar/Polyester.  Therefore,  laminates  having  the 
same  ply  orientations  as  the  panels  were  tested.  Appendix  A 
details  how  the  material  properties  were  determined.  The 
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properties  used  for  the  analysis  of  Kevlar/Pol yester  panels 


are 


E1  =  6.7932  *  106  psi 
E2  =  0.44502  *  106  psi 
t'12  =  0.50369 
G12  =  0.1782  *  106  psi 

Graphite/Epoxy  and  Kevl ar/Pol yester  cylindrical  panels 
were  analyzed  using  STAGSC-1  to  predict  their  collapse  load 
under  compression.  The  dimensions  of  the  panels  used  in  the 
analysis  are  illustrated  in  Figure  2.1.  The  top  curved  edge 
was  clamped  and  allowed  only  a  displacement  in  the  x 
direction.  The  bottom  curved  edge  was  fully  clamped.  The 
vertical  edges  were  simply  supported,  allowing  rotations  and 
restraining  any  ci rcumferential  and  radial  displacements 
along  the  boundary.  Mathematically,  these  boundary 
conditions  may  be  expressed  as 

Top  edge:  u=free;  v=w=w, x=w, y=0 

Bottom  edge:  u=v=w=w,  =w,  =0 

x  y 

Vertical  edges:  u=w,  =free;  v=w=w,  =0 

x  y 

Once  the  boundary  conditions  are  set,  the  next  thing  to 
consider  is  the  finite  element  mesh  size  to  be  used  to 
insure  a  correct  answer  is  converged  upon.  Becker  [2] 
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conducted  a  study  to  determine  the  optimum  mesh  size  for  the 
finite  elements  of  composite  cylindrical  panels.  He  found 
that  1/2"  x  1/2“  elements  were  ideal  for  the  panel 
geometries  under  consideration  using  the  QUAF  410  or  QUAF 
411  elements.  Also,  Wilder  [1]  comments  that  accurate 
results  are  obtained  with  STAGSC-1,  using  QUAF  elements,  if 
the  element  size  is  between  1/2"  to  2/3”  square.  As  a 
result,  1/2“  x  1/2"  elements  were  used  for  all  analysis 
performed  in  this  thesis. 

The  last  thing  to  consider  before  the  analysis  can 
begin  is  the  type  of  element  to  be  used.  The  QUAF  410 
element  was  chosen  for  use  in  the  analysis  done  in  this 
thesis.  It  was  felt  that  since  the  vertical  edges  were 
supported,  the  moderate  rotations  encountered  could  be 
adequately  modeled  by  the  410  element.  The  411  element 
would  be  the  element  of  choice  for  large  rotations.  Another 
reason  for  selecting  the  410  element  is  the  difference  in 
the  resulting  computational  time  between  the  use  of  the  two 
elements.  An  analysis  using  the  411  element  with  its  eight 
additional  degrees  of  freedom  can  take  up  to  four  times 
longer  than  an  analysis  using  the  410  element. 

Table  4.1  shows  the  collapse  load  for  each  material  and 
ply  orientation  studied  assuming  a  geometrically  perfect 
panel  with  the  above  stated  boundary  conditions.  The 
experimental  loads  that  were  shown  in  Table  2.3  ranged  from 
24-62%  below  these  analytical  loads.  These  large 


discrepancies  can  partially  be  explained  by  considering 
geometric  imperfections  of  the  curved  panels. 


Table  4.1 

STAGSC-1  Collapse  Load  of  Undamaged  Panels 


[0/90] £S 

[±45/90/0]q 

[0/903 3. 

Load 

(lbs) 

Gr/Ep 

5926 

4652 

10124 

Kev/Pol y 

21 1  7 

1504 

5524 

STAGSC-1  Imperfection  Analysis 

The  manufacture  of  geometrically  perfect  cylindrical 
panels  is  nearly  impossible  to  achieve.  Janisse  [3]  in  his 
work  showed  that  imperfections  do  exist  in  panels  and  can 
have  a  profound  effect  on  their  load  carrying  capability  for 
certain  ply  orientations.  The  imperfections  he  measured 
were  on  the  order  of  the  thickness  of  a  single  ply  (0.005"). 

STAGS  has  a  built  in  capability  to  impose  a  geometric 
imperfection  into  the  panel.  The  imperfection  is  constant 
through  the  thickness  of  the  panel  and  of  the  form  [22] 

WQ  =  WAMP  *  cosj"(x-X1)  *  k/Xl!  *  cos["(y-Y1)  *  k/YlI  (4.3) 
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where 

WAHP  =  Amplitude  of  Imperfection 

XI, Y1  =  Location  of  Center  of  Imperfection 

XL , YL  =  Half-wavelength  of  Imperfection  in 
x  and  y  direction  respectively 

The  question  that  must  be  answered  at  this  point  is, 

"What  imperfection  pattern  will  best  model  the  actual 

imperfect  panel?”.  To  determine  the  appropriate  pattern, 

two  ply  orientations  were  examined,  [0/90]2s  and 

[±45/90/0]  .  The  material  analyzed  was  Graphite/Epoxy.  It 
s 

is  assumed  that  the  Kevl ar/Pol yester  panels  would  respond 
similar  to  the  Graphite/Epoxy  panels  since  the  imperfections 
are  geometrical,  not  material  related.  First,  consider 
[0/90] 2s  Graphite/Epoxy  panels.  Table  4.2  shows  the 
collapse  loads  calculated  by  STAGSC-1  from  a  nonlinear 
collapse  analysis  of  panels  with  imperfection  of  the  type 
described  by  Equation  (4.3).  Four  different  imperfection 
patterns  were  analyzed. 

(1)  0.005"  with  one  half-wave  in  each  direction. 

(2)  0.04"  with  one  half-wave  in  each  direction. 

(3)  0.04"  with  three  half-waves  in  each  direction. 

(4)  0.005"  with  three  half-waves  in  each  direction. 


4-10 


Table  4.2 

Imperfection  Analysis  of  [0/90]2s  Gr/Ep  Panels 
(Collapse  Load  of  “Perfect"  Panel  =  5926  lbs) 


Imperfection  Amplitude 


#  of 

half  waves 

0.005" 

o 

o 

Load 

1 

6334 

4995 

(lbs) 

3 

2880 

6516 

The  surprising  result  of  this  analysis  is  that  for  two 
of  the  imperfection  patterns  studied,  the  panel  actually 
strengthens  as  a  result  of  the  imperfection.  To  understand 
this,  consider  the  shape  of  the  "perfect"  panel  at  the 
collapse  load  as  calculated  by  STAGSC-1  shown  in  Figure  4.2 
Figure  4.2  actually  shows  the  deformed  shape  of  the 
"perfect"  panel’s  centerlines. 

First  consider  imperfection  pattern  (1).  The  geometry 
of  this  imperfect  panel  is  very  close  to  that  of  the 
“perfect"  panel.  Janisse  [3]  showed  in  his  work  that  the 
imperfect  panels  tended  to  displace  in  the  same  pattern  as 
the  imperfection.  If  one  thinks  of  the  imperfection  as 
being  imposed  on  the  panel,  then  for  this  first  case,  we 
have  a  panel  that  wants  to  deform  into  3  half  waves 
longitudinally  and  5  half  waves  circumferentially  being 
forced  to  deform  into  one  half  wave  in  each  direction.  One 
can  see  that  it  will  require  some  extra  energy  to  pop  the 


(b) 


Shape  of  "Perfect 
(a)  Longitudinal 


Figure  4.2 

[0/90] Gr/Ep  Panel  at  Collapse 
Shape  (b)  Ci rcumferential  Shape 
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extra  waves  through  until  the  panel  deforms  into  one  half 
wave  in  each  direction.  Hence,  the  collapse  load  is  higher 
than  that  of  the  perfect  panel. 

Next,  consider  imperfection  pattern  (2).  Here  we  see 
that  as  the  amplitude  of  the  imperfection  is  increased,  the 
panel  weakens  somewhat.  With  an  imperfection  as  large  as 
0.04",  the  applied  load  not  only  causes  membrane  stresses 
but  also  begins  to  impart  some  bending  to  the  panel.  This 
additional  bending  weakens  the  panel. 

Third,  consider  imperfection  pattern  (3).  Again  we  see 
the  collapse  load  increases  beyond  that  of  the  "perfect" 
panel.  In  this  instance,  the  imperfections  are  large  enough 
and  the  curvatures  in  the  panel  change  often  enough  that  the 
panel  responds  similar  to  a  folded  plate.  The  effective 
longitudinal  length  of  the  panel  is  decreased,  causing  an 
increase  in  the  collapse  load. 

Lastly,  consider  imperfection  pattern  (4).  The 
collapse  load  of  2880  lbs  is  within  3*  of  the  experimental 
collapse  load.  This  imperfection  pattern  is  closest  to  the 
deformed  shape  of  the  "perfect"  panel  determined  by 
analysis.  Therefore,  less  energy  is  required  to  cause  the 
panel  to  collapse  because  the  imperfections  present  in  the 
panel  are  in  the  direction  the  panel  wants  to  naturally 
move.  The  exact  geometry  of  the  imperfections  is  not  known 
and  in  general  is  random  and  can  vary  from  panel  to  panel. 
However,  the  pattern  of  three  half  waves  longitudinally  and 


3 


circumferentially  seems  to  model  the  actual  imperfections 
very  well  by  smearing  them  out  evenly  over  the  entire  panel 
geometry. 

Consider  now  the  same  four  imperfection  patterns 
imposed  on  a  [±45/90/0]  Graphite/Epoxy  panel.  Table  4.3 

S 

shows  the  collapse  loads  computed  by  STAGSC-1  for  these 
cases. 


Table  4.3 

Imperfection  Analysis  of  [±45/90/0]  Gr/Ep  Panels 

S 

[Collapse  Load  of  "Perfect"  Panel  =  4652  lbs) 


Imperfection  Amplitude 

#  of 

half  waves 

0.005" 

o 

o 

Load 

1 

4692 

4930 

(lbs) 

3 

4065 

6200 

We  see  the  same  trends  here  for  the  [±45/90/0] 

9 

orientation  as  for  the  [0/90]2s  orientation  with  one 

exception.  For  imperfection  pattern  (2),  STAGSC-1  computes 

a  greater  collapse  load  than  that  of  a  "perfect"  panel. 

Recall  the  collapse  load  was  lower  for  the  [0/90]„ 

zs 

orientation.  The  two  ply  orientations  see  the  load  applied 

in  the  same  way,  but  the  [±45/90/0]  orientation  has  a 

s 

greater  bending  stiffness  than  does  the  [0/90]2g 
orientation.  The  ±45°  plies  are  the  reason  for  the 
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increased  bending  stiffness.  This  also  accounts  for  the 
dramatic  increase  in  the  collapse  load  of  a  panel  with 
imperfection  pattern  (3).  The  angled  plies  help  resist  the 
bending  moments  in  addition  to  the  geometry  decreasing  the 
effective  length  of  the  panel. 

One  last  aspect  of  the  imperfection  analysis  to  comment 
on  involves  the  sensitivity  of  the  panels  to  initial 
geometric  imperfections.  We  see  by  the  analysis  that  the 
ply  orientation  determines  the  sensitivity  of  the  panel  to 
imperfections.  [±45/90/0]  was  least  sensitive  to 

S 

imperfections  while  [0/90]„„  was  most  sensitive.  The  author 
feels  that  the  ply  orientation  is  critical  in  characterizing 
a  cylindrical  panel’s  imperfection  sensitivity. 

Now  that  the  appropriate  imperfection  pattern  has  been 
determined,  we  perform  the  analysis  for  each  ply  orientation 
and  material.  An  imperfection  was  located  at  the  center  of 
each  panel  with  three  half  waves  in  each  direction.  The 
amplitude  was  set  equal  to  0.005".  Thus,  the  following 
parameters  are  substituted  into  Equation  (4.3). 

XI  =  6.0" 

Y1  =  19.0986° 

XL  =  4.0" 

YL  =  12.7324° 
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Figure  4.3  shows  the  shape  of  the  imperfection  along  the 
center  of  the  panel  in  the  x-di rection.  Table  4.4  shows  the 
collapse  loads  of  the  panels  accounting  for  the 
imperfection. 


Table  4.4 

Collapse  of  Imperfect  Panels 


[0/90] 2S 

[±45/90/0]s 

[0/9°]3s 

Load 

Gr/Ep 

2880 

4065 

7473 

(lbs) 

Kev/Pol 

833 

1320 

2338 

While  the  imperfections  that  are  actually  in  the  panel  are 
not  exactly  modeled  by  this  sinusoidal  imperfection,  the 
experimental  loads  are  within  3-16X  of  the  loads  in  Table 
4.4. 

Figures  4.4  -  4.9  shows  the  Load-Displacement  curves 
for  each  panel.  The  displacements  are  in  the  direction  of 
the  applied  load.  STAGSC-1  curves  for  a  perfect  and 
imperfect  panel  are  shown  as  well  as  experimental  data. 

The  curves  show  the  same  features  for  each  case 
studied.  Note  that  the  STAGSC-1  curves  for  a  perfect  and 
imperfect  panel  show  some  post-buckl i ng  strength  for  the 
[0/90] 2S  and  the  [0/90]3s  ply  orientations.  However,  the 
experimentation  was  halted  after  the  panel  buckled  as  it  was 
not  the  intention  of  the  research  to  investigate 
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Load  (lbs) 


Load  (lbs) 


5000.00 


Figure  4.5 

Load-Displacement  Curves  for  [±45/90/01  Gr/Ep  Panels 


Load  (lbs) 


2500.00 


Figure  4.7 

Load-Displacement  Curves  for  [0/9012s  Kev/Poly  Panels 
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Figure  4.8 

Load-Displacement  Curves  for  [±45/90/0]  Kev/Poly  Panels 
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Figure  4.9 

Load-Displacement  Curves  for  [0/90l3s  Kev/Poly  Panels 


post-buckling  strength  of  cylindrical  panels.  We  see  that 
the  imposed  imperfection  pattern  has  little  or  no  effect  on 
the  initial  stiffness  of  the  structure.  In  fact,  no 
difference  in  stiffness  is  noticed  until  the  collapse  load 
of  the  imperfect  panel  is  reached.  One  other  thing  that  we 
did  not  expect  to  see  was  the  large  difference  in  the 
stiffness  between  the  two  analytical  panels  and  the  actual 
panel . 

Analytical  vs.  Experimental  Stiffness 

Looking  at  Figures  4. 4-4. 9,  one  notices  the  slopes  of 
the  analytical  curve  are  much  steeper  than  that  of  the 
experimental  curve.  The  finite  element  model  is  thus 
stiffer  than  the  actual  panel,  which  is  expected.  However, 
the  amount  of  the  difference  is  somewhat  surprising. 

On  first  glance,  one  must  suspect  the  model  being  used. 
Since  boundary  conditions  play  a  role  in  the  stiffness  or 
flexibility  of  a  structure,  changing  the  boundary  conditions 
was  considered.  A  STAGSC-1  analysis  was  performed  on 
[0/90] 2S  and  [±45/90/0]g  Graphite/Epoxy  panels  allowing  the 
ci rcumferential  displacement,  v,  to  be  free.  Figures  4.10 
and  4.11  compare  the  analytical  imperfect  panel,  with  both  v 
fixed  and  v  free,  and  the  experimental  panel  for  the  ply 
orientations  considered. 

As  expected,  changing  the  boundary  conditions  changed 
the  response  of  the  panel.  However,  little  or  no  change  is 


apparent  in  the  initial  stiffness  of  the  structure  when 
considering  the  in-plane  displacement,  u.  For  the  [0/90]2g 
panel,  no  change  can  be  seen,  where  as  for  the  [±45/90/0] 

S 

panel,  a  slight  increase  in  flexibility  is  noticed.  Figures 
4.10  and  4.11  show  that  while  boundary  conditions  do  have  an 
effect  on  the  stiffness,  it  is  still  most  likely  the 
imperfections  in  the  panel  which  cause  the  increased 
flexibility.  Even  though  the  two  analytical  curves  have 
imperfections  modeled  in  them,  they  do  not,  by  any  means, 
model  the  actual  imperfections  exactly.  Therefore,  the 
author  feels  it  is  the  actual  imperfections  of  the  panels 
that  cause  the  increased  flexibility. 
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5.  CONCLUSIONS 


Based  on  the  analysis  and  experimentation  performed  in 
this  thesis,  the  following  conclusions  may  be  drawn. 

1.  For  composite  cylindrical  panels,  there  is  a 
minimum  impact  energy  below  which  no  damage  is  sustained  by 
the  panel.  This  minimum  energy  depends  on  the  ply 
orientation  and  the  thickness.  The  quasi-isotropic 
orientations  are  the  most  easily  damaged  where  the  thicker 
cross-ply  orientations  require  more  energy  to  cause  damage. 

2.  Photographs  and  C-Scans  show  that  the  damage  area 
appears  to  remain  relatively  constant  for  a  range  of  impact 
energies.  However,  we  can  not  say  the  total  damage  remains 
constant  because  we  don’t  know  how  the  damage  propagates 
through  the  thickness  as  the  impact  energy  increases. 

3.  Panels  impacted  with  energies  between  1-3  ft- lbs 
show  little  decrease  in  compressive  strength  (2-6*)  as  a 
result  of  the  impact.  However,  the  delamination  caused  by 
the  impact  spread  as  a  result  of  the  compression  test. 

Thus,  impact  damage  can  be  quite  catastrophic,  even  though 
the  damage  may  not  be  visible,  particularly  if  the  structure 
is  subjected  to  fatigue. 

4.  Simulating  impact  damage  by  implanting  teflon 
inserts  as  Wilder  [1]  did  may  be  used  to  get  an  idea  of  the 
effect  impact  has  on  the  global  collapse  load  of  the  panel. 
Wilder’s  technique  would  yield  conservative  results  since 
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the  teflon  inserts  provide  a  totally  debonded  region.  The 
damaged  area  may  still  have  some  strength  left  in  it 
depending  on  the  amount  of  energy  imparted  to  the  panel. 

This  technique  has  only  been  tested  on  quasi-isotropic 
laminates.  This  conclusion  can  not  be  extended  to  cross-ply 
or  any  other  laminates  without  further  testing. 

5.  Imperfections  must  be  accounted  for  when 
determining  the  collapse  load  of  undamaged  cylindrical 
panels.  This  is  especially  true  for  cross-ply  laminates. 

Not  accounting  for  imperfections  led  to  differences  in 
analytical  vs.  experimental  loads  of  24-62%. 

6.  The  imperfection  pattern  of  3  sinusoidal  half  waves 
in  each  direction  (longitudinal  and  ci rcumferenti al ) ,  with 
an  amplitude  equal  to  the  thickness  of  one  ply  (0.005"), 
seems  to  most  accurately  model  the  random  imperfection 
present  in  a  8"  x  12“  cylindrical  panel. 

7.  Quasi-isotropic  ply  orientations  are  relatively 
insensitive  in  the  longitudinal  direction  (x)  to  geometric 
imperfections  where  as  cross-ply  laminates  are  extremely 
sensitive  to  imperfections. 

8.  Changing  the  boundary  conditions  on  the 
circumferential  displacement  (v-fixed  vs.  v-free)  had  little 
effect  on  the  flexibility  of  the  panels.  Thus,  it  is 
primarily  the  actual  imperfections  of  the  panel  that  cause 
the  larger  than  expected  displacements  in  the  u-direction. 
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9.  Microphotographs  of  the  cross-section  of  damaged 
composite  panels  show  a  great  deal  of  cracking  and 
delamination  through  the  thickness  of  the  laminate  not 
visible  from  C-Scans  or  stereo  X-rays.  Information  gained 
from  microphotographs  will  be  essential  to  develop  a  model 
to  predict  damage. 

10.  The  assumption  that  E./E  ,  E0/E  ,  and  i>  /i>  are 
the  same  for  Kevlar/Pol yester  as  they  are  for  Graphite/Epoxy 
was  valid.  The  resulting  equivalent  orthotropic  material 
properties  were  within  5 %  of  the  measured  properties  of  the 
woven  laminates. 
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Appendix  A :  Material  Properties  of  Kevl ar/Pol vester 


Kevlar/Polyester  is  a  weave  material  that  cannot  be 

made  into  unidirectional  coupons  like  Graphite/Epoxy. 

Therefore,  determining  the  material  properties  is  not 

straight  forward.  Also,  STAGSC-1  can  not  analyze  a  weave 

material,  so  an  equivalent  orthotropic  laminate  had  to  be 

developed  to  model  the  Kevlar/Polyester. 

Kevlar/Pol yester  coupons  with  ply  orientations  of 

[0/90]2s  and  [0/90]3s  were  tested  in  tension.  The  results 

from  testing  these  coupons  yielded  E  ,  E  ,  and  i>  for 

x  y  xy 

Kevlar/Polyester.  G12  was  determined  from  the  f±45l2s 
coupons.  From  the  tensile  tests, 

E  =  3.6371  *  106  psi 

A 

Ey  =  3.6371  *  106  psi 

t'  =  0.061936 
xy 

G12  =  0.1782  *  106  psi 

Next,  in  order  to  calculate  E1  ,  E2,  and  i>12  for 
Kevlar/Polyester,  a  relation  between  these  laminate 
properties  had  to  be  assumed.  The  assumption  is  that  the 
ratio  of  laminate  properties  to  unidirectional  properties 
for  Kevlar/Pol yester  is  the  same  as  for  Graphite/Epoxy. 
Therefore  the  effective  laminate  engineering  properties  for 
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Graphite/Epoxy  were  calculated  for  ply  orientations  of 


[0/90]2s  and  [0/90]3g  using  the  following  equations  [23]. 


Ex 

=  (AnA22 

-  A12)/hA22 

Ey 

=  (a11a22 

-  A*2)/hA,, 

(A.l  ) 

l>xy 

=  A12/A22 

where  the  A.  . ’s  are  elements  of  the  extensions!  stiffness 

*  J 

matrix  and  h  is  the  laminate  thickness.  Equations  (A.l)  are 
only  valid  for  balanced  ply  orientations  (equal  numbers  of 
+0  and  -6  layers)  symmetric  about  the  middle  plane.  For 
[0/90] 2s  and  [0/90]3s  Graphite/Epoxy, 

Ex  =  10.9549  *  106  psi 
Ey  =  10.9549  *  106  psi 
i’xy  =  0.0385 

In  order  to  obtain  an  initial  guess  for  the  equivalent 
orthotropic  engineering  properties  of  Kevlar/Polyester,  we 
now  use  the  assumption  that  the  ratios  E1/E)<,  E2/Ey,  and 
t'l 2/i'xy  were  the  same  for  both  materials.  Using  these 
ratios,  the  following  equation  was  used  to  calculate  E1  for 
Kevlar/Pol yester . 


where 


0  = 


J  Gr/Ep 


(A. 3) 


Similar  expression  were  used  to  calculate  E2  and  t>12.  The 
values  obtained  were 


E1  =  6.7932  *  106  psi 
E2  =  0.44502  *  106  psi 
t>1 2  =  0.50369 
G12  =  0.1782  *  106  psi 

Using  these  values,  E  ,  E  ,  and  v  for 

x  y  xy 

Kevlar/Polyester  were  calculated  using  Equations  (A.1)  and 
compared  to  the  experimentally  determined  laminate 
properties  as  a  double  check.  The  calculated  values  are 
within  5%  of  the  experimental  values.  Thus  the  above  values 
for  E1 ,  E2,  i'12,  and  G12  were  used  in  STAGSC-1  for 
Kevlar/Pol yester . 
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Appendix  B:  Classical  Laminated  Plate  Theon 


Since  the  present  study  deals  with  composite  panels,  a 
brief  review  of  classical  laminated  plate  theory  is 
presented.  Interested  readers  are  referred  to  Jones’  [24] 
text  for  a  detailed  presentation  of  classical  laminated 
plate  theory. 

For  an  orthotropic  material  in  a  state  of  plane  stress, 
the  stress-strain  relations  may  be  written  as 
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The  Q.’s  are  the  reduced  stiffnesses  and  are  functions  of 
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the  material  constants. 
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Equations  (B.1)  are  defined  in  the  principal  material 
directions.  However,  these  directions  often  do  not  coincide 
with  the  body  or  laminate  axis  system  natural  to  the  given 
structure.  Therefore,  a  coordinate  transformation  must  be 
made.  Transforming  from  the  material  axis  system  (1-2 


B- 


directions)  to  the  structural  axis  system  (x-y  directions) 
we  obtain 


Q  Q  Q  1  f  e 

_lt  _12  _16  x 

Q  Q  Q  <  e 

_12  _22  _26  y 

Q  Q  Q  Y 

16  26  66  J  l  ' xy 

where 

Q4.  =  Q„«cos*0  +  2(Q,_  +  2QC,  )s i n20cos20  +  GLosin*0 

®i2  =  ^Qn+  Q22-  4Q66)s,in20cc,s20  +  Q12(sin*0  +  cos40) 

Q22  =  Qltsin40  +  2(Q12  +  2Q66 )sin20cos20  +  Q22cos40 

516  =  (Qn“  q12-  2Q66) si n0cos30  +  (012-  Q22+  2Q66 )sin30cos0 
°26  =  (Q:i‘  Qi2‘  2Q66)sin30cos0  +  (Q12-  Q22+  2Qg6 )si n0cos30 
°66  =  (Q11+  °22’  2Q\z~  2Q66)sin20cos20  +  Q66(sin40  +  cos40) 

( B.  4 

The  angle  0  is  the  angle  between  the  1-axls  and  the  x-axis 
shown  in  Figure  B.l. 

Knowing  how  the  stresses  and  strains  vary  through  the 
thickness  of  a  laminate  is  necessary  in  order  to  define  the 
extensional  and  bending  stiffnesses  of  the  laminate.  Here, 
we  assume  the  individual  laminae  are  perfectly  bonded 
together,  infinitesimally  thin,  and  are  non-shear 
deformable.  The  displacements  through  the  thickness  are 
continuous  so  that  no  laminae  can  slip  relative  to  another. 
We  also  make  the  following  assumptions,  collectively  known 
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Figure  B.1 

(a)  Fiber  (1-2)  Axis  System  (b)  Transformation  into 
Structural  Axis  System 


as  the  Kirchhoff  hypothesis  for  plates  or  Ki rchhoff-Love 
hypothesis  for  shells. 


(a)  A  line  originally  normal  to  the  middle  surface  of 
the  laminate  remains  normal  after  the  laminate  is  deformed 
(i.e.  shear  strains  7  and  7  are  negligible). 

xz  yz 

(b)  Normals  have  constant  length  (ez  =  0). 


The  strains  for  the  laminate,  as  long  as  they  are 
assumed  to  be  small,  may  be  written  as 
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where  the  middle  surface  strains  are 
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and  the  middle  surface  curvatures  are 
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In  Equations  (B.6)  and  (B.7),  u  ,  v  ,  and  w  are  the 

OO  O 

displacements  of  the  middle  surface  in  the  x,  y,  and  z 
directions  respectively. 
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Substituting  Equation  (B.5)  into  Equation  (B.3)  we 
obtain  the  stresses  in  the  kth  layer  of  the  laminate  in 
terms  of  middle  surface  strains  and  curvatures. 


f  y 

f 

r  o  i 

r  y 

o 

e 

K 

X 

r  -i 

X 

X 

0 

IO 

II 

4 

e° 

•  +  Z  < 

K 

► 

y 

L 

l 

y 

X 

r 

K 

v  sy  J 

k 

l  xy  J 

l  xy  J 

J 

( B .  8 ) 


The  resultant  forces  and  moments  acting  on  the  laminate  are 
obtained  by  performing  the  following  integrations  through 
the  thickness  of  the  laminate. 
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N  is  the  number  of  layers  in  the  laminate  and  zfc  and  zkl 
are  shown  in  Figure  B.2.  Figure  B.3  illustrates  the  forces 
and  moments  acting  on  the  laminate. 

Upon  substitution  for  the  stress  vector  of  the  kth 
layer  and  performing  the  integrations,  we  obtain  for  the 
forces  and  moments 
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The  A,  B,  and  D  matrices  are  symmetric  and  are  defined  by 
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